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In this paper, we study hydrodynamics and thermal behaviors of micro-heat exchanger of different cross 
sections, with hydraulic diameters (Da) of 0.4—0.8 mm and Reynolds numbers (Re) of 300-900. Validity of the 
model is proved by comparing simulation results based on classical Navier-Stokes equations with experimental 
data and Tuckerman correlation. Effects of dynamic viscosity, hydraulic diameter, Reynolds number and cross 
section of the micro-channels are investigated. The results indicate that the variable viscosity effect should be 
taken into account, and rod bundle micro-channels have the highest heat transfer coefficient of the three types 


of micro-channels. 
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I. INTRODUCTION 


Compact heat exchangers based on micro-channels show 
a nice foreground over miniature nuclear power system, and 
learning flow behaviors of fluid in micro-channels is of im- 
portance for designing micro-fluidic systems. 

Investigations [1-5] have been devoted to flow character- 
istics of water in micro-channels. Koo and Kleinstreuer [6] 
studied viscous dissipation effects on convection through cir- 
cular and rectangular micro-channels and concluded that vis- 
cous heating was a function of the hydraulic diameter, the 
channel aspect ratio, and the Brinkman number, and that ig- 
noring viscous dissipation could affect accurate flow simula- 
tion in micro-channels. 

Judy et al. [7] experimented with distilled water, methanol, 
and isopropanol in round and square micro-channels fabri- 
cated from fused silica and stainless steel, and measured fric- 
tional pressure drop in hydraulic diameters of Dy = 15- 
150 um over Reynolds number of Re = 8-2300. The results 
showed that the values of f Re agreed well with the conven- 
tional Stokes flow theory over the tested range of Reynolds 
numbers. 

Gamrat et al. [8] performed both two- and three- 
dimensional numerical simulation to analyze thermal en- 
trance effects on micro-channel convection. They found that 
the continuum model of conventional mass, Navier-Stokes 
and energy equations were adequately accurate to simulate 
the fluid flow and heat transfer for the micro-channels. 

Against researches in the past decade in a micro-channel 
of single shape, however, there is just a handful of studies, 
experimental or theoretical, on flow characteristics of micro- 
channels in different shapes. In this paper, we study hydro- 
dynamic and thermal behavior of fluid flow in micro-channel 
heat exchanger of different shapes. The effects of dynamic 
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viscosity, hydraulic diameter, Reynolds number and cross 
sections are investigated. 


Il. EXPERIMENTAL SETUP 


The experimental system used in this study is divided into 
three parts—the test section, the water driving system, and the 
dynamic data acquisition section, as shown schematically in 
Fig. 1. 


Data acquisition 


Reservoir 


Collection 
container 


Bypass valve 


P: pressure sensor T: thermocouples 


Fig. 1. (Color online) Schematic of the experimental system. 


Effective hermetically sealing and heat insulating measures 
were taken for the test section. Fluid flow through the test sec- 
tion could be adjusted by the control valve. Deionized water 
was used as the working fluid. A flow meter was connected 
at the entrance of the test section. Two pressure sensors were 
fixed on the inlet and outlet ducts to measure pressure. Type- 
K thermocouples were selected to measure temperatures of 
solid wall and liquid. Heating pillars powered by a DC power 
supply, which can be regulated, were used to supply constant 
heat flux. The measurement signals were transmitted to a 
computer. 


The physical configurations of the micro-channels heat ex- 
changer are shown schematically in Fig. 2. The dimensions 
of cooling passages are given in Table 1. 
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TABLE 1. Dimensions of the micro-channels (um) 
Dn A WwW Ry Re W2 W: L 


400 2200 220 200 210 140 20 10000 
600 3300 330 300 300 200 20 10000 
800 4400 440 400 420 280 20 10000 


II. MATHEMATICAL MODEL 


According to the literature, Navier-Stokes equations are 
still valid for laminar flow of water in micro-channels of hy- 
draulic diameter as small as 0.4mm. The model is derived 
from continuum based on conservation equations of mass, 
momentum and energy, with the following basic assumptions: 
(1) negligible effect of gravity and other forms of body forces, 
(2) incompressible Newtonian fluid, and (3) steady laminar 
flow. 

The governing equations are Eqs. (1)- (3), the equations of 
continuity, momentum and energy, respectively. 


VV =0, (1) 
pV- VV) =-Vp+ uV? V 4S, VT -defV, © 


prep,s(V VT) =k VT + 4, (3) 


where, py is water density, p is pressure, V is the velocity, 
kf is the water thermal conductivity, ks is the silicon thermal 
conductivity, p is the viscosity of the water, c, ¢ is specific 
heat of water, T is temperature. 

For the solid, Eq. (3) can be written as: 


ksV°T =0, (4) 
with viscosity-temperature sensitivity as: 
Su = du/dT, (5) 
and the viscous dissipation term as: 
¢ = pdefV /(udefV /2). (6) 


The viscosity is calculated by the parameters supplied by 
the NIST in the temperature range of liquid water, given as: 


u = 0.11157 — 9.51523 x 10-47 + 2.7249 x 10-67? 
— 2.61107 x 107°T®. 


(7) 
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Fig. 2. Schematic diagram of computational domain. 
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The boundary conditions at the inlet, the outlet, the heating 
surface and the fluid-solid interface are given by Eqs. (8)— 
(11), respectively. 


u = Un, V = 0,w = 0,T = 293K, (8) 
P; = Pyy(latm),v = 0,w = 0, (9) 
— ks0T/Oy = qon, (10) 
Tar = Týr, 
uļp = 0, vip = 0, wļp = 0. 


The governing equations with boundary conditions were 
solved by the commercial CFD code, FLUENT. The equa- 
tions were discretized by means of a fully implicit second 
order finite volume method with modified upwind advec- 
tion scheme. The SIMPLE algorithm is used to resolve the 
pressure-velocity coupling. The numerical solution process 
is shown in Fig. 3. 


Built governing equations 
Built initial and boundary conditions 
Mesh and generate computing nodes 


Establish discrete equations 


Discrete initial and boundary 
ondition 


Given control parameters 


Solve the discrete equations 


onvergence 
or not 


Display and output results 


Fig. 3. The solution process. 


IV. EXPERIMENTAL RESULTS 


The experimental and simulation results for the f Re and 
heat transfer coefficient as a function of the Re number, are 
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Fig. 4. Comparison of the experimental and simulation results for f Re (a), and heat transfer coefficient as a function of Re (b). 


plotted in Fig. 4. The trends of experimental data and pre- 
dicted curve are in good agreement over the entire range of 
Reynolds number being studied. 


To compare with experimental data, the local thermal re- 
sistance is calculated using 


R(x) = (Tmax (£) — Tin)/4, (12) 
where, R(x) and Tmax(£) are thermal resistance and maxi- 
mum temperature at x cm from the entrance, Tin is inlet water 


temperature, and q is heat flux at the heating area. 

Validation of the numerical method was performed under 
uniform heat flux condition over the two side faces, as a 
benchmark case. Fig. 6 shows the thermal resistance data 
measured in this work and results by Tuckerman [9]. It can 
be seen that they agree well. The slight differences may 
be caused by different physical properties and outlet condi- 
tions. Apparently, the numerical method is appropriate for 
the present study. 


Thermal Resistance (cm?:K:W-1) 
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Distance from the entrance (cm) 


Fig. 5. Comparison of the measurement results (0) of thermal resis- 
tance with the data (O) by Tuckerman [9]. 
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Fig. 6. The local f Re (a) and local Nusselt numbers (b) as a function 
of LT. 


040601-3 


WU Ge-Ping, WANG Jun and LU Ping 


mY -Æ Rectangular 
0.8 -¥= Rob bundles 
7 i -^ Round 
O vr 
© 06 
Ge 
S 
D 0.4 
'> 
LL 
0.2 
0.0 


0.4 0.5 0.6 0.7 0.8 
Hydraulic diameter (mm) 


Fig. 7. Friction factor vs. hydraulic diameter for micro-channel of 
different shapes. 
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Fig. 8. Effects of heat flux uniformity on the f Re (a) and thermal 
efficiency of micro-channel of different shapes (b). 
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V. DISCUSSION 


As shown in Fig. 6(a), the local f Re we simulated is in 
good agreement with analytical solution by Yilmaz [10]: 


ea 13.76 | 64) + K/L* — 13.76/V L+ 
€e)Yilmaz = VLF “Í + 0.98 x 10-4 K3-14/(L+)?’ 
(13) 
L* = L/(ReD,). (14) 


The f Re reaches a constant value of about 120 when the 
hydrodynamic flow develops fully with a constant viscosity. 
However, the viscosity varies actually due to increasing water 
temperatures, and the f Re decreases continuously, as shown 
in Fig. 6(a). Therefore, the viscosity effect should be taken 
into account. 

The calculated local Nusselt numbers agrees well with the 
theoretical results by Shah [11], too, as shown in Fig. 6(b). 
The data of variable viscosity are always less than the data of 
constant viscosity, as the viscous dissipation decreases heat 
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Fig. 9. Thermal efficiency vs. D, for three types of micro-channels. 
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Fig. 10. Nusselt number as a function of Da under uniform or non- 
uniform heating condition. 
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Fig. 11. Effect of Reynolds number (a) and hydraulic diameter (b) on heat transfer for the three types of micro-channels. 


transfer efficiency. Differences in two sets of data disappear 
gradually with the flow development, as the viscous dissipa- 
tion effects decrease. 

The friction factors of laminar flows, as a function of hy- 
draulic diameter, are shown in Fig. 7 for micro-channels in 
rectangular, round and rob bundle cross-sections. For the 
same cross-section, the friction factor decreases with increas- 
ing Dp. The rob bundle micro-channels have the highest 
friction factor, which causes the greatest energy dissipation. 
On the other hand, for D, = 0.6-0.8 mm, the round micro- 
channels have the smallest friction factor. 

Figure 8(a) shows that the f Re rises with Reynolds num- 
bers under both uniform or non-uniform heat flux, but the 
f Re curve of non-uniform heating is always higher than that 
of the uniform heating, and the differences increase gradually 
with Re. 

The performance index, defined as the ratio of heat trans- 
ferred to total pumping power, is introduced for better under- 
standing of the overall performance of the micro-channels, 


Q = PrCp (Tout) = Ti 
Ap : 


Nett = (15) 


Ppump 

Figure 8(b) shows the performance index, as a function of 
the Reynolds number, of differently shaped micro-channels of 
the same hydraulic diameter. Generally, the performance in- 
dex decreases with increasing Reynolds number. This is due 
to the increased pressure drop which results in higher pump- 
ing power as the flow rate increases. 

The performance index of different types of micro- 
channels with the same Reynolds numbers is illustrated in 
Fig. 9. The performance index decreases with increasing hy- 
draulic diameter. Of all types of micro-channels, the round 
type has the highest performance index. The rod bundles type 
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